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ABSTR?VCT 


Diiaocyanates# as well as some triols, which are 
reactants for polyurethane systems, are known to exhibit 
unequal reactivity. Polyurethane network formation has been 
modelled in this work taking into account unequal reactivity 
due to tho asymmetry and induced asymmetry of the functional 
groups in the monomers. Using the model developed, stoichio- 
metric ratios of reactants leading to gelation in polyurethane 
systems have been calculated by employing a combination of 
kinetics of roactlons and probabilistic methods. The effect 
of secondary reactions of allophanate and biuret formation 
was also evaluated. It was shown that unequal reactivity can 
explain the occurrence of a minimum in geltime below an 
isocyanate to hydroxyl ratio of 1 , 

Unequal reactivity models wore also used to compute 
CO2 ovolution and average molecular weight for polyurethane 
foam systems. The agreement between observed and computed 
rates of CO2 evolution was satisfactory qualitatively. 



CHAPTER 1 


INTRODUCTION 


Polyurethanes arc callaJ the polymers of the 80* s. 

Itiis is beccaiso they have a wide range of useful properties 
and they hove not been fully exploited yet. Presently new 
uses aro being found and the annual consumption of this group 
of polymers is growing very fast. A look at the wide variety 
of uses of those polymers will explain this* 

Mg id ■ Polyuroth ane Fo am s are used for construction, 
packaging, insulation, automobile bumpers end some special 
uses, 

2* Floxlblo Foams are used for cushions, packaging and some 
special uses, 

3. Polyurethane Elastomers are used in extra strong automobilo 
tyros, automobilo fascia and shoe soles. Ihe automobile 
industry found this polymer very attractive and analysts 
predict that in future the door panels and marry of the car 
body parts are going to be made of polyurethanes* 

The reasons for such usefulness are many 

1) Mochanical properties of the polyurethanes can be varied 
from soft olastomors to rigid character by simple alterations 
in the initial composition of reactants* 

2) Manufacturing time of the polyurethanes can be made very 
small as the reactions are very fast compared to other polymers 
even at ordinary temperatures, inf act by using some tin 



cotalysts tho .jolation timo £t»r the polymur can bo brought 
diwn to as lew as a fuw socc'nds. 

3) Other Pix 'duct properties such as dimensional stability, flamo 
rutardation, oxidation resistance, color stability and other 
desirable properties can bo imparted to polyurethanes by 
addition of a number of suitable additives, 

4) L^rMe .-ner-r/ savin is in prc'ducing polyurethanes, can bo 

of foe tod by Reaction Injection Molding, This new technique of 
mcilding usos tho fast nature of the polyurethane formation 
reaction. Such fast nature allows pumping of less viscous 
monomers into the mould where polymerization takes place in situ. 
Thus cost of molting and pumping of viscous polymer is avoided. 
The savings are considerable as calculations show^ that petroleum 
oqulvolant a feedstock plus energy required for production is 
tho least for Reaction injection Molded pcljfurethane among all 
the synthetic matters including aluminium, steel and other 
mineral products* 

Polyurethanes^ are so named becouse the linking 
reaction between the reactive groups of the moncxners ( Hydroxyl 
and isocyanate) leads to the formation of urethane bonds as 
given below j 

O H 
il 

-v^OH + O e C = N -> —'0 _ C N 

( isocyanate) ( urethane) 

Usually bifunctional isocyanates are used, the most common 
ones being mdI (4,4* -Diphenylmethane-diisocyanate) , 

TDI ( Ibluene diisocyanate) and HDI ( Hex am ethylene diisocyanate). 



The hydroxyl compounds can be a mixture of diols and triols. 

The diols are usually polyothers or polyesters terminated with 

OH groups and are referred to as polyots. Besides the above 

laain reaction^ many side reactions (including that with foaming 

agents) also occur. The final properties of the crosslinked 

polymor depend upon the extents to which all the reactions 

took place and how they contributed to the formation of the gel, 

3 4 

Literature survey revealed some reaction models and simulation 
but none which accounts for all the simultaneously occurring 
reactions and their complex nature. 

While it has been reported that diisocyanates exhibit 

very complex unegual reactivity^ only simple unegual reactivity 
5 

models have been used in published works to model this system. , 
Moreover the complexities arising duo to the presence of 
reactive blowing agents have not been accounted for. The aim 
of the present work is to account for these complexities and to 
compare theoretical predictions with experimental data whesre I 
available. i 

j 

The models presented in this work would be used in ; 

Chapter 2 to predict the compositions leading to gel formation, 

I 

I 

Crosslink density, urethane group concentration and solfraction | 
equations can be found in Appendix C, such predictions are 
useful in calculating the mechanical properties of the | 

polyurethanes, 

I 

In many applications, the gelled polymer is taken ; 

6 

out of the mold and further cured using secondary reactions i 



4 


( allophanate and biuret formation) . The effect of these 
side reactions is studied in Chapter 3, 

The models have been used to predict gel times end 
were checked with available experimental data. Those results 
are contained in Chapter 4. 

For the production of foams it is necessary that 
the GO 2 evolves at the right time when the polyurethane 
has gained sufficient viscosity. The viscosity is dependent on 
the average molecular weight. Hence it is necessary to look 
at the CO 2 evolution rate as well as the molecular weight of 
the reaction mixture as a function of time. These calculations 
are contained in Chapter 5, 



CHAPTER 2 


A typical network forming polyurethane foam system 

B 

has as initial reactants, a tricl (B— ‘ — ^B) , a diol (B-B), a 

diisocyanato (A-A) and water as blowing agent (C-C), From among 

these the isocyanates and water exhibit unagual reactivity. 

Three basic types of unequal reactivity have been identified. 

7 8 

structural asymmetry or asymmetry refers to the situation 
whore the difforait functional groups which are chemically 
similar react at different rates. Those cases where the 
reactivities of all the functional groups of a monomer are equal 
but different from that of the very same functional groups 

present at the end of a polymer chain have been termed induced 

8 9 

asymmetry or first shell substitution effects • Cyclic 

8 

monomers which gaierate a chemically equivalent group after 

reaction constitute the last category of unequal reactivity. 

It is possible that a monomer may simultaneously possess one or 

more of the basic unequal reactivity characteristics, several 

isocyanate monomers show the simultaneous presence of structural 

2 lo 11 

asymmetry and first shell substitution effect ' * , while 

hydroxyl containing compounds may exhibit structural asymmetry* 

BASIC REACTIONS ; 

In polyurethane systems the reaction between hydroxyl 
group and isocyanate group produces the urethane linkage: 

H O 

-w^N = C s= 0 4- HD'-^ — > ~~N - C - 0--^ (1) 
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Tho reactions between water and isocyanate group are: 

H O 

i 

N = C = 0 + H 2 O ~~N-C-OH ( 2a) 

H O 

— - N - ^ - OH — ^NH2 + CO 2 ( 2b) 

The gaseous CO 2 so liberated in the second reaction acts as 

the foaming agent. The amine generated in the second reaction 

can react with isocyanate group producing polyurea linkages: 

H O H 
i H ^ 

M = C = O + H 2 N — ' — N - C - N ( 3) 

Apart from these, allophonato and biuret formation reactions 

can contribute to the chain branching. While some catalyst 

systems may promote these side reactions, generally their 

2 6 

effect can be ignored except at high temperatures ' , In 

the present chapter these side reactions will be assumed to be 

absent. 

BASIC MODEL 
Urethane Formation ; 

Let the di isocyanate monomers be denoted as A^A 2 thereby 
explicitly recognizing that the two monomeric isocyanate groups 
will react at different rates. Let the diol be denoted by B 2 and 
the polyol by B^, The same symbol B is being 



7 


assi'.. Oil wo ?■’ 1 hydroxyl groups as all of them have been assumed 
to bo equally reactive. The structural asymmetry of the 
diisocyanato mononer can be represented by the following 


equations : 


Vi + 


*■ 1^2 ® 



ic 

ii B 
^ 2 


(4) 

(5) 


The and present at the und of a polymeric chain have been 

clearly distinguished ty putting a star as they will react 

at a rate different from that of the corresponding monomeric 

group duo to first shell substitution effects. The first shell 

substitution effect can be modelled as 

* 

^1 

- - ( 6 ) 
k* 

- (7) 


Art A' 


2^1 


B 


Aj_A2 4- 


B 


Linear polymerisation of such A^A 2 moncmer with only B 2 monomer 
has been analysed by Gendhi and Babu^^ and their results can 
be applied to analyse urethane formation by reactions (4) to 
(7). 

Reactions with water ; Since reaction (2) is a very fast 
one# it will bo assumed to occur instantaneously. Then a water 
molecule reacts with one isocyanate group and# in effect# 
generates the amine group directly. The amine garoup can then 
react with one isocyanate group ty reaction (3), Such a 
sequence of reactions is similar to those of a cyclic monemor 
and in particular those between an anhydride and hydroxyl group. 



0 


In other words water can be modelled as a cyclic monomer. Let 
water be denoted by CD. The C group reacts first consuming 
one isocyanate group and generates a D group (amino group). 

The D group does not exist in the monomer and hence can not 
react till the "adjoining" C group reacts. Therefore the 
D group can be present only at the end of a chain while the 
C group cannotform the end of a chain. The D group can then 
react with one isocyanate group. The C and D groups can react 
with the isocyanate groups at different rates. For simplicity it 
will be assumed that the asymmetric characteristics of A^A 2 
monomer are inherent to the ^■^^2 the same character- 

istics used in urethane formation are applicable here also. 

In other words we are assuming that the relative reactivities 
of polymeric and monomeric and 1^2 groups are unaltered 
irrespective of the group with which they react though their 
absolute reactivities may be different. The reactions with 

water may thoi be written as follows: 

Kk, 


* 


'^2^1 

+ CD 

^ a^a^cd 

(8) 



Kk , 

(9) 

^ 1^2 

+ CD 



The A^ and A 2 groups of the products of reactions (8) and (9) 


are polymeric in character. Thus we will write 

* 

* Kk^ 


A2A^ + CD 


* 


Kk, 


A.A„ + CD 


The D ( amino) group reactions are 






• AaA^CD 


A^A 2^D 


A^A^ D- 


( 10 ) 

( 11 ) 


( 12 ) 



+ D-' 




^2 * 


'k 




;-'x 

^R^ ^2 




^1^2 ° "' 


V^l D 


■'-'•Ai^.2 D 


(13) 

(14) 

(15) 


Reactions (4) to (15) then form the basic model of the system. 


PROBAEilLITY OF NET^^^ORIC FORMATION 

Initially let only mononers ^2-' ®f ' present. 

For the system being considered^ gelation occurs through formation 
of structures like 

- B •'■ 

- B B - A^A^CD. A^A^ 

.. Q " 

If tho probability of occurrence of the above structure regard- 
less of the number, of intervening bifunctional CD and BB units, 

i€ 

or the branching coefficient , can be calculated, the gel 

l/ 

point can be predicted folla-i/ing Plory , The probabilities of 
finding structures of tho type 

BA2^A2""/ '■*~‘A2^A2CD’"' ■, cA^A^BB— ^ i~A^A2B— 

are needed to calculate tho branching coefficient. All such 
probabilities have boon show^ to be related to the 

fraction of the total number of the various groups that have 
been consumed in direct reaction. The probabilities and their 
interrelations can be calculated once tto.e kinetics of the 
reactions is known. An example of such aoalculation is given 
in Appendix 



This can bu 


Probability o£ finding unrcactcd Aj^A 2 monomer t 
simply dcfinju as tho ratio of unreacted rnonjmer concentrations 
at any given time and initially. 


(X, = 


A / Va' 

U -o 


Probabilities of reaction of A groups ; During 'the course of 

reaction moncmeric groups arc lost by reaction (4) as v?ell as 

by ( 5 ^ , Reaction (4) consumes a monomeric group by direct 

reaction whereas reaction (5) converts a moncmeric group to 

Let j- 

a polymeric A 2 group, ^ A^ be the moles of moncmeric A^ groups 

that have been consumed by reaction (/^-i^^pcr unit volume at any 

I’o 

time. Following Gandhi and Babu let us define the probability 

that a monomeric A^ group has n-acted as 

r 


a: 


= 


\ A *' + A 


(16) 


Such a definition has tho advantage that it readies a value of 
unity at the end of the reaction when A^?k 2 monomer is in 
stoiciiicmetric deficiency. Similarly the probability that a 


moncmeric ^2 Stroup has reacted is 


a: 


''2 - 


fA’ + a; 


Note that 


o = ■>- t 


r 


A^ A, 


(17 ) 


(18)' 


When A 2 moles of moncmeric A^ groups have reacted, the same 

* . 
number of polymeric or A^ groups have been created, IF LA^.j 

is the ccncontration of unreacted polymeric A^ groups at 

y* ^ 

any time, then A 2 - !. A^ - number of polymeric A^ groups have 



react.' r'^ection (6), Hence, the probability o£ reaction o£ 


polymeric grcups can be defined as 




A- — ■, A, i 
/ - 1 - 

A^ 


(19) 


Similarly, the probability of reaction of polymeric A^ groups 
is defined as 

.r 

— = (20) 


Ai 




Probabilities of reaction of B groups ; Let ■ ®2--o 

represent the initial concentrations of the polyol and diol 
respectively. Let \ B<beth€ concentration of unreacted B ends 
present at any time. The probability that a randomly selected 
B group has reacted is defined as 


B \ 


f + 2 


(21) 


2-0 


Prob^llltl<s. of of . CD units : Let fC'.^ and 

respectively represent the concentrations of the CD monomer 

initially and at any time. Let D represent the concentration 

of unreacted D ends at any time. The number of CD units that 

have reacted per unit volume is C » C , Hence the 

o 

probability of finding a CD unit where C group has reacted can 
be defined as 


■C ' - -.c ' 

N' = — (22) 

•PJo 

During the course of reaction, (C<^ - 'iC j moles of CD monomer 
units have reacted and hence that many moles of D ends have been 
created per unit volume, out of theseonly \ D"\ moles per unit 



volumo of unrcactod D ends o.r.-e present, Honoe the moles per 
unit volume of CD units in which both c end D groups have reacted 
is eciual to ‘ c' - .c • - 'dj. The probability of finding a CD 
unit vjhe.re both C and D have ruactod can bo defined as 

'C «• 0 ^ H) 

S = — (23) 

' C 

o 

Again note that ^ 'f. and •' will attain unity at the end of 

the reaction if monomer is in stoidaicxnetric excess. 


Stoichiometric rolcitionship i The total nirmber A groups that 
have reacted must be equal to the total number of B, C and D 
groups that have reacted. Hence 

S = (1 -tcA*) + aJ (1 o + ^ 5 

In terms of the above defined probabilities^ the first equality 
reduces to 

S s 'A :Ki (1 - ^2 ^ ^ ^ ^ + -a ^ ^ - -^1 ) (1 + (x )•’ ' Aq /Ml )'( 1 - A ^ )'j; 

(24) 

Let us define the fraction of poly functional groups/^^ , 




£ I B,'- 
- f o 


£ + 2rB„;, + 2 >C. 


(25) 


■ f ' o " 2^' o ■ o 

and fraction of hydroxyl groups contributed to the diol, d^ 

2 i'3, ■ 


d = 


2’o 


^ >£^0 + 


(26) 


and the fraction of blowing agent groups, b, 

2 -C ’i 


D 


^i®f’o + 2 -CM 


b 


(27) 



The above three quantities are related through 


r = (1 o- d) (1 - b) or d(l - b) = 1 - - d (28) 

be defined 

Furtlior let ' • the stcichicfnetric ratio, R,/as the ratio 

of initial concentrations of A groups and the ccncentrations of 
all the groups with v7hida A can react 


R = 

Then using (25 ) to 


2 (29) 

f.Bjio + 2 :b2K_, + 

(29), tlie total B,c,D groups that have 


reacted, s, ic 


given by 


R(1 c d) 


%+ (1 


)(l+9 




(30) 


Probability of a structure of the type«-*A^A 2 BBc- t The structure 
being considered can bo formed in four ways j 


(i) Monomeric group reacts, then the polymeric A 2 group 

reacts with a B 2 unit and the B end reacts subsequently, 

★ 

can represent such formation (ii) Monaneric A 2 group 

reacts, the polymeric group reacts with a B 2 unit and the B 

* 

end reacts aftonv-ards, SB- roprxisents this type of 

formation (iii) One B group of a B 2 unit can react, the other 
B end can then react v7ith moncmcric group and the polymeric 

ic 

A 2 group can react thereafter, **^2^1 symbolises this type of 
formation (road fran right to left), (iv) A B 2 unit can react, 
monomeric A 2 reacts with the other B end and then polymeric A^^ 

ic 

group reacts, -A^A^BB- is the representation of this type of 
formation. The total probability of formation of -A^A^BB- type of 
unit, irrespective of the internal arrangements of A^A 2 or the 
ways in vftiich the unit was formed, is the sum of the probabilities 



of ocaarronce of the four typws of sequences of events listed 


dijove*. 

The probability of occurronGCj of the first of the above 
sequences will be calculated as an cxanple. The prob^ility 
is given by 


(probability of finding) 
( a monomeric reacted ) 

(group next to a bond ) 


(probability that) (probability that) 
(polymeric A2 reacts) (the B end has ) 
( with a BB unit ) i^e^-ctcd 


•.r 

(-^) ( 


2 ; Bo’ .t’ 
. Z.' O'"' 


) ( ) 


Similarly the probability of the second, third and fourth 
sequences can be calculated to be, respectively. 



2iB 


2 '.o 



•iB 2 ‘ 



,* ^'^ 2 :. oh 

'2' s 




■* 


S ' *1 


Hence the total probability of finding -A^A^BB- unit Pgg is 


given by 
Pbb « 2 





s 



■* 

■''1 


2 ’.B 


2 ** o 


-.2 


After substitution of ( 24 ) and ( 30 ) 

i< ic 

^ *0 -’^1 ^ 


Pbb - 4 


^l >2 


2 '“'! 


■<1(1 



Z'-'o+d-d-'OCf'+i) 

I • 

(31 ) 


Probability of a structure of the type - Aj^A2CD- : Analysing 
in a manner similar to the previous section, there are eight 
ways in which the structure being considered can fom. Using a 
notation similar to the earlier section these eight ways can be 

^ V . 'At 'A A 

listed as : -A^A^CD- , -A^^A^DG- , -A2^l'2D- , -A^A^DC- , Aj^A2C]>- , 

^ t • < A A 

— A^A2DQ- # — A^A^^CD— f — A^A^DC- * As an example, the fifth 
structure may be interpreted to have formed when a monomeric A2 



grojp roacfcs x^?ith CD moncmur and subsequently the polymeric 


group erio iroup gunorated i. jact. Followiri^, the* method 


employed in the earlier section, the probabilities of the above 

aF 




.-.ight structures can bo calculated and are given by ° 


'^1 ,* * C ^§ 

— =— o<. --> of. 2 

S * ^2 * ^1 ~ ' 






s ^ 

^ . _^2 QS' 

s ^ s 


c ^ 

C "** 


•i 0<.*2 . 


^ A2^ * 

~ , — 1 respectively. The total probability o 

S ^ 


n'F 


finding -A^A^CD- structure irrespective of the internal arrange- 
ments of the A^A 2 and CD units is then the sum of all the 

above eight probabilities or 


PcD 


Y* 

Ar + Ao c ' 

1 z 2 1 o 


In terms of probabilities and other definitions adopted earlier 

■k 


^CD ~ 


- 1 2 - -'2^ 

(1-''2)(1 +* 2 ) +C2(l-^l)(l+ 3L^ 


(1 - d -r) C 


2'^t+ ( 1— '■ ( X + S ) 

(32) 


Probability of a structure of the type -7vj^A2P-.' : There are four 

ways in which this structure can form: 


The probabilities of these four structures respectively are 

.r 


a: 


-i 

S • 2 


* f.B.! :b 

* f O ' 


2 


* * '®£Pl='' 


S • 1 


f Tb 


f- O 


- r 
"2 


* 

• " 1 * 


fiB.. 

f- o 




1 - 

^ , ■*'■ 2 * The total 


prob^ility of finding this structure Pg is given by 

■“f 

S * S 


and after substitution of earlier definitions 



4 



o<3_Tr^2)(i+^ )+<^(i -o<^ ) (i+i>5_ ; • 9p^a--d^f) ( V +S’) 

( 33 ) 


Br':indijLnq Coefficient ; Consider the following structure 


The probability of the ctoove structure is simply given by 

However the -A^7>.^CD- and -iv^A^BE- units can be interchanged and 
it will load to (m + n)l/(mini) number of stanjctures with 
the same probability. Hence the total probability of finding 
such a structure irrespective of the internal rearrangements is 


(m + n ) t 

mi nj 



<''CD> 


m 




Branching coefficient is equal to the probability that a B 
group selected at random frem one of the B^ units will lead to a 
B group of another B^ unit rcgardlessof the number of ccxinecting 
bifunctional units* Thus branching coefficient is given by 

CO CO 



(m + n )1 

ml nl 


1 - - 


p 

"■BB 




or 


13 


Following arguments of Flory , the gel point occurs when 


(f-1) 


1 _ p „ p 

CD BB 




( 34 ) 


The above criteria along with kinetics of the reactions can be 
used to predict the fractional conversions at which gelation 



will occur, “Hac inter-relations between probabilities have 
been in Appendix I assuming secend order kinetics. The 

gel criteria thv..n simplifies to 


(f-1) 


4 fi (1 - d) 


^ l-d-y=) (f +S J-4 ) ^1-d-p J 

(35) 




The above criteria must be simultaneously satisficKi along witii 
the stoichicmetric relation which now reduces to 


0(1 


1 ^* + k,x<iK)]= r( 


1+K^C<1 

Both these equations are proved in the appendix I* 


j^[2f^+(l-d-/)Cf+$ )] 


(36) 


RESOLTS AND DISCUSSION 

The gelation criteria given by (34) can be used to 

predict whether a particular initial composition specified by 

R, f-' and d will lead to the formation of an infinite network or 

not. Further if a gel wore to form, the fractional decrease in 

the initial functional group concentrations required before 

gelation cai occur can also bo predicted by (34), ResCilts have 

chapter 

been obtained in this e,>ar^for reaction of a diisocyanate vjith 
a triol (f = 3 ), a diol and water. 

Recovery of appropriate limits : Under appropriate restricticais, 
the present model must reduce to those developed by earlier 
workers, A few such examples are given to demonstrate the 
validity of the present model. The kinetics presented in the 



<i.ppGndix I will bo usod as aiid whon nocossary 


(i) Limit of oqual rtJ activity ; When ^ 2^1 — ~ 

it ^ 

^I'^l = ^ ^^ 2/^2 — ^2 ” both the functional groups, 

, 3 ncl P^ 2 t equally reactive in their monaneric as well as 
polymeric form. Further if K = 2 and = 0.5, the CD monomer 
behaves like Dp montmaer , Flory's equal reactivity limit 
must bo rocoV'.„rcd with these restrictions. The following 
relationships can bo easily established from the appendix with 
the above restrictions ; 


ir= 





0^1 


l+o<3_ 


-n 






1 


(37) 

If is defined as the fraction of total A groups that have 
reacted, thcan 


Pa 




a: 


:(i-H3<^) 


2 



Q<i 

l+0<]_ 


ic 

(1 +'■>^ 2 ^ 


From (37) and (38), it can bo shown that 


(38) 


0(1 P ;^(2 - p^) ^ 

_= 


Pa 


*Pa 


Similarly if is defined as the fraction of B, c and D groups 
that have reacted, then, it can be shown that 

* 'So <T-^S) 

^i?£'o + ^'®2Jo ^ (?]o 


Pb 


= 


(39) 


If all the above results are substituted in the stoichiometric 
^relationship (36), it reduces to R = pg and the gelation 
criteria (35 ) reduces to 



.V 1 


(40) 


( £- 1 ) 


■>' Pa % 

( l--( I- j ) P; 5 j?Q 


Both these results arc identical to that given by Flory^^. 

Note that Flory's r is 1/R and BB monomer in Flory's reference is 
the <Vk monomer of this paper. 


( ii) Limit of only first shell substitution effects This 
limit is obtained when = 1, To make CD equivalent 


to 62 # we let K = 2 and = 0,5, Hence 


.2 


. 2 


■1 - ' 2 ' '"I 


■'2' 
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Miller and Macosko considered gelation of 
Ag with B 2 where exhibits first shell substitution effect. 
To compere therefore we have to put g = 2 in their results and 

I - Q 

f = 2 , k = 1 , d = O, b = O in ours. Miller and Macosko define 
to bu the fraction of the initi^A^ units in which i number 
of A groups have reacted. It is easy to show that 


p = ^ = i— 1 

° 1+ k- 


' O i-t- A I 

r .* ^ 


Pi 


iiv 

r A ’ 

o 


2 1- 


1 + 


'1 


and 


2 A 


1 2 


•' 1"'2 


rK 


1 + './■ ■ 


If those results are substituted 


So * 

in the gel criteria given by Miller ard Macosko, it reduces to 
★ 

2 / (1 + ■^■' 2 ) = 1 * For the case being considered^ the 

2 

stoichiometric relation (36) simplifies to +^>^. 2 )/(l + ^ 1 )=^ 


and the gel criteria (35) is given by 2 l + -">^ 12 ^ = 

these two are combined the gel criteria is identical to that 

given by Miller and Macosko f 

rj. 

Bokarc and Gandhi considered the crosslinking of f 
functional epoxide with g /2 functional primary amine, where the 
amine hydrogens exhibit first shell substitution effect. They 
also considered the side reaction of the epoxide group with the 



s--condary hydroxyl group gonoratt-id. If the hydroxyl group is 

considorod to bo nonroactivc and if g = 2 , their results aro 

fi^r the reaction between with v\hor 0 only first 

shell substitution effect. As hydroxyl groups are not reactive 

only secondary amino S-j_ and tertiary amino T 2 (defined in 

12 - 

refcrence 9) exist and gel criteria given by them' reduces to 
2(f-l) ~ "^2 number of tertiary amine units or 

in our notation, A 2 _A 2 units in which both and A 2 have reacted. 


Hence, 


T 2 2 



_ p ^2 

1 -i- 


and gel criteria given in equivalent to 2( f-1) 1 + = 1, 

This is identical to that predicted by (35) and (36). Note 
d = b = O since the reaction is between ^]_A 2 and only. 


(iii) Limit of only structural asymmetry; This limit is 

i‘c ic 

obtained when =: K 2 = 1, K = 2, = 0,5, Hence - ^ . , 

.2 2 5 
\ =s 2.' - " , b = ; • For this case Miller and Macosko have 

given the gelation criteria to be r = %2 • 1 + (f -2)0^ : . 

and %2 fraction of A^ and A 2 groups consumed and are 

easily, shown to bo given by ( 1-f /( 1 + 

'•^1 t^A +'■’^ 2 ^'^ ^ respectively. Noting that r = 1/R and 

a£ != /-, the gelation criteria of Miller and Maaosko reduces to 

^ tf-2)/] /(I -f K^9<.3_) = 1, Simplifying ( 35) 

and (36) it can bci shown that (35) is identical with that of 

Miller and Macosko, 


Gelation Envelope : For a given system, the reaction rate constants 
are fixed and (35) can predict the fractional conversions at 



which gelation will occur onco R, , ‘ and d arc spocified. 

Hov;^'\rir, the uiaxintum possible conversion is ccmploto conversion 
of those graaps which are in stoichiometric deficiency. For 
example When R < 1, gel can form only when all.x/. 1, Hence for 
a given ' and by gelation condition .and stoichiometry, 
there exists a n minimiom value of stoichiometric ratio, 

at which gelation can occur when all = 1, For any R '' 
gelation is not possible for the same value of r' and d because, 
before the required amounts of poly functional B groups can 
react ^^1^2 would be exhausted. Similarly there 

will be a such that for R " Fj^^x^ gelation will not occur. 

To illustrate tliis point consider as an example gelation of 
with and B 2 , The system is one of equal reactivity and 
d = { 1 - ' ) since b = O, The gel criteria has been shown to be 


1 


^ Pa Pb 
-( 1 -; ) 



(40) 


and the stoichiometric equation is given by Rp^ = Pb* R C 1, 

the maximum value p^ can attain is 1 while pg is always less 
than 1, Using stoichiometry, the godation criteria siit^lifies 


to 


2 R p 


A 


1 - ( 1- i- ) 


— 2 

Rp,t 


■"y 


or R 


"T 

Pa 


(1 + F ) 


Thus when p^j^ s= 1, R takes its minimum possible value? 

^in “ 1 +. 1 "* ^ • Clearly gelation can not occur at the same 

if R <" Rj^jLn that would require ' 1. If ^/\iin 9<^l3‘‘tion 

can occur but for P;^ <C. 1- Similarly if R > 1, the maximum 
value possible for pg is 1 while Now the gelation criteria 

will reduce to 



1 


Equal reactivity 

%/ ’ 


t: < 

CL :■ 

: a/>!nttrM-i nr\f nriccShfo 



FIGURE 1 } Gelation envelop© for the case of equal 
reactivity. (K»2, Kjj,=Kj^bK2«1) 



2 j 

(R - (1 





or R ^- ( 1 + p ) 



Hen CO tho maximijm value of the stoichiometric ratio, at 

which gelation can occur is given by (1 + f ) when has been 
oquatoc?. to one. Gels will form only if R 4. ^\ia 3 { 

^ ^ point. Thus gelation is possible 

only if ^ ■^. E^ax ^ given r , Those two limits then 

define the boundaries of tiio gelation envelope and it is shewn in 
figure 1 for equal reactivity case. Gelation willoccur only 
if the initial corttposition falls inside the envelope. The 
j"' and vs, ^ will be referred to as the lower and upper 

boundaries of the envelope respectively. The gelation envelopes 
for different parametric values of the reaction rate constants 
and d will bo presented in this suction. First to be considered is 
the case when d = O for simplicity, i 


Gelation envelope in the absemoe of bifunctional mononner; 


( i) EJ^in '^^^sus ! This gives the lower 

boundary of the gelation envelope, would be less than unity 

and therefore monomer is in stoichiometric deficiency. 


The value of is obtained when all A groups have been 

consumed or when = i. -ihen (35) reduces to 


4 0^/' 


2j-b +( u-pii-l,) 


^ 1 


or 


ti 4 P 
J 1 r 


2Pe+ (f-S )1^ ) 


(41) 


and (36) reduces to 

= 2ff6 + (1 -p ( '{ + S ) 


2R 


(42) 



In tho above we have taken f = 3 and d = O. By definition 
is always positive and . Then (41) can be rearranged 

to get 

(4- ^ _ 2£: ):^{1 - .-O ( ) 

and hence the minimum value that £: must attain before gelation 
con cccur is given by 

^ - 2 r. O or : j- ■> 0.5 

Further note that gelation will occur at = 0,5 if :“ = 1 

I" 

regardless of the values of the reaction rate constants. 

The value of ut J'= 1 therefore is obtained from (42) to 

be 0,5 and this result is canmen to all gelation envelopes. 

If gelation wore to occur at . 3 = 1 and since, 3 = 1 a.lso 

moans that ^, = x = 1^ the only value possible for R . is 

min 

unity ty (42), Conversely vhen R_j_ = 1/ £ = = 1 and ciny 

mxin i'*' 

value of 3 ,,. 0 will load to gelation by (41), But when 
;■ = O only = 1 will satisfy (41) and (42), Therefore, 

R^in = 1 when o forms a common point for all gelation 

envelopes • 

Since is completely consumed, only the 

time rotjuired for gelation and not the lower boundary of 
the gelation envelope is affected by the reaction rate constants^ 
ratios and K 2 , Only the reactivity ratios K and 

which describe the effect of the behaviour of the blowing agent, 
influence the results. 

Effect of K; This parameter indicates the effect of 

the relative reactivity of the c group. Let 

. be considered 

the two obvious limits given by K = o and K->go^ When K = 0, 

the blowing agent is nonreactive and hence acts sinply as a 



diluont. Both'/” and ^ arc thercf:.rG soro and (41) predicts 

that gelation will occur as seen as ^ rs^aches 0.5 and hence by 

(42), result indicates that as J-’ is 

increased, more /Vi mcn.omor rfxjuirid for gelation. It is an 

expected result since as is increased, the fraction of the 

rti^octivo mononer increases wiiile the fraction of nonreactive 

blowing agent dccrt.ases. As K — > eg the blowing agent CD 

reacts instontan'-ously and completely before any B groups can 

be consusiicd. Hence 'Y' = S = 1 snd therefore gelation will 

again C'ccur as soon as = 0,5, By (42) therefore 

^\nin — ^ j^instly this result shows that for the same 

value of requiirod for K = O is less than that for iC-— . oo. 

This is to be expected since CD does not consume any A 2 A 2 

when K s O while CD requires stoichicffnotrioally equivalent of 

Aj^A 2 viien K-— ^co. Secondly the above result ^ows that, 

unlike when K = 0, as is increased, A^A 2 required for gelation 

decreases. This can be understoodsince CD mononer coisumes 

A^A 2 monomer at the beginning itself ferrming VJiCD- type prepolymer 

but not resulting in crosslinking. Hence as P is increased, 

the CD monomer fraction decreases and thus the requiromont of 

A 2 ^A 2 is reduced. 

The blowing agent can therefore act in two opposing 
ways: as a nonreactive diluent and as a chain extending consumer 

of AA, As is increased, the former effect tends to increase 

while latter effect tends to decrease R . , The values of 

mm 

/ 

K and J-' determine the effect that dominates. It is pertinent ' 

here to point out that, had CD been assumed to be equally reactive 

as the theory would predict a raonotcnic dec 3 a:ase in R . as P is 
^ min J 



increased. Numerical results illustrating the effect of K are 


presented in figure 2 for K.j^ 1. For K > 2, the diluent 

>.;ffect is n..-t there at ail and mnnotonically decreases 

v.'i. t)i inero3.sir-'-; ;■ , However V7hcn K is suffieiontly low^. such 

as 0,1^ very interesting behaviour that reflects a balance of 

ttie two -.if foots alixiady described is observed, At such a low 

value of K, the blowing agent reacts very slowly and should 

b.-have liK._ a ''iluont. Hut at lew valuos of ■ , largo amrunts 

of CD arc present and thus the low reactivity is compensated 

i-r by the higli concentration, As a result, a significant amount 

vf AA is consumed in chain extension, H^nc^-j at low values of , 

the blowing agent behaves like a chain extending agent and 

incroasod values of cause a decrease in R . , However at 

min 

sufficiently high values of , large; ameunts of fast reacting B 

grouyjs arc present to consume AA me nomer leaving the blowing 

agvjnt o'ssontially unr^acted , The blowing ag-.jnt now essentially 

plays the rolv: of a diluent and hence increased results in an 

increased This conclusion is confirms by calculations 

which show that ij = O in curve 2 till the minima is reached. 

To tho left of the left boundary of this curve, gelation will not 

take place since too few molecules are there to create 

sufficient number of connected brandh points. No gelation is 

possible to the right of tho right boundary since excess of 

B^ mcnomor is present leading to too many unconnected B groups. 

All the curves of figure 2 indicates, as expected, that 

at constant J-' an increased value of K would mean more consumption 

of AA by CD or larger chain lengths between br.ancii points and 

that can be 

hence increased R Another conclusion ./ - drawnfrem these 






cur'/'-Gis that, i£ R is fixed, gels can be formed even if large 
amc'unts if blowing agont arc added if a retarding agent for 
it:; reaction is £-.un>‘. 


Effect :;.f Again, it will be fruitful to investigate 

the two 1 iinits of = Q and ■ ct . When K„ = O, the D 

K iv 

group is unrcactivo and hence CD acts like a mono functional 
blocking agent and =0, The gel criteria now reduces to 


r 




'to in 

!Z 


(43) 


2 f 


and stoidhiornetric relation is given by 


2R = 2 f ;S + ( 1 p ) Y' ( 44 ) 

When /3 = 1/ i is also' equal to unity and R = 1 + 0,5 by 
(44), We can therefore infer from (43) that gelation will 
occur only i£ p '11/3, This is shewn in figure 3 in curve 1 • 

The* top portion of curve 1 was obtained for c<<l and i 

fo'r the indi exited values of the parameters. The relations (43) 
and (44) can not bo combined into one analytical opression 
and numerical results can be presented only after assuming a 
value for K, The portion of the curve 1 for 1/3 ^ been 

calculated for K = 2, As is increased, the fraction of 
blowing agent is decreased and the amount of AA monomer 
wastcfully consumed and blocked decreases. Thus decreases 

as / is increased as seen in figure 3. 

w* 

When ^ CQ the D grioup reacts instantaneously or in 
effect each c greup that reacts consumes two A groups. Thus 
'Y = and therefore gelation occurs as soon as = 0,5. 

Stoichiemetry r equires that 





R =: 


Ml-; )■:= 1/. )a -d-p 


-sK'l 


-in:l tliorcfcro 


K 


R 


min 


4* ( 1 — H) j 1 ■!• ( — 


Results cf ilhis liiait aro reprosentod by curve 4 o£ figure 3 
for K = 2. Note that at p =0,^= 0,5 is a meaningless 

criteria and hence gel will form only if this 

limit also, like in the limit of K--^ o. as the value of is 
increased^ ‘\nin and is to be expecto.l as the blowing 

agent behaves like a chain extending agent. 

The effect of increasing K„ on R at constant ^ is 

R man 

quite unexpected as increases at first and decreases later* 

This can be observed from figure 3* At low values of the 
D group does consume A groups Ujt at such slow rates that CD 
monanor essentially behaves like a blocking agent. Hence, 
increasing K^, from zero effectively leads to increased consumption 
of A groups but not to effective chain extension as long as 
is sufficiently low. Thus initial increases in from zero 
will increase However as is increased further the 

blocking effect vanishes and CD will new participate in 
effective chain extension or connecting branch points. This 
effect will then decrease the requirements of A groups and hence 
^'min decrease. Another useful conclusion from these 

results is that if is increased to sufficiently high values 
by catalysis, gels can still be obtained by addition of larger 
amounts of blowing agent at constant R, 

Figures 4 and 5 display the combined effects of K and 
for different sets of values. 







( ii) versus - : This gives the upper boundary 

e£ tlio gelation envelope, ^ax greater than unity or 

the A gr-.ups are in stoichiometric excess and hence the upper 

boundary is obtained \fthen ■' = / = ;>= 1, The only parameters of 

significance tli- r-fore are K_ , K- and Ko as (35) and (36) 

ik* JL Z 

now simplify to 



1 + K 


A 


"^2 + 



(45) 


ani 


-^1 


1 + K, 


1 ^2 ^ ^ +d'v. 




R 


(46) 


■'1 


In particular note that d does not have any effect. The only ^ 
point common to all gelation envelopes is p = 0 and R^^x ~ 

This can bo soon from (45) since if p = O, (45) can be satisfiet 

T*f ik 

only =-'.2 = 3. ■v\hich in turn implies that = •a 2 = 1# 

MumuriCcil results illustrating the effect of and 


are shown in figure 6* 

i 

Effect of Note that K^p/ 1 since labelling of 1 

and 2 on A groups is arbitrary. For simplicity let us first 

consider the case when = 1, If the value of is 

greater than unity^ A-, reacts faster than A and hence if largo 

1 

amounts of A^A^ monemor are present, all B groups are consumed 
by A 2 groups leaving groups unrfsacted. To promote crossl inkin 
the A^ groups have to be f orced to react and this can be achieved 
only by decreasing R at constant p , Thus the effect of increase 
value of is to decrease as shown in figure 6 by curve 2. 


Effect of IC^ and K 2 on The 

values of the above parameters for the 
various curves are: -1: 1,1,1,; 2:5, 1,1 
3:1, 4,4; 4:10, 10,10 and 5: 2.5, 10,10 



Sff ;..ct of or K^; First consider the case when = 1 
and = K2. If cr K2 is greater than unity the polymeric 
A groups react faster than monomeric I\ groups or the processes 
of connection of branch points and chain extension are promoted, 
Hcncelarger amounts of AA monomer can be utilised or is 

increased at constant . This is shown by curve 3 of figure 6 , 
At con'-tant thcrofuro the effects of increased K, 
and ( or K2) on in the opposite directions* 

Thu net ^-ffect is determined by the numerical magnitudes of these 
constants. This is illustrated in figure 4 by the other curves. 
It is observed in all these curves that as is increased 
^.max When '• is increased the fraction of 

increases. Gelation will not occur if monomer terminates 

all 3 ,C/D groups without linking two branch points and this 
happens when A2^A2 is in excess of However the chances 

of two branch points being ccrmectod arc increased when the 
fraction of the poly functional groups or r-' is increased and 
thureforu mere A2_A2 monauer can bo tolerated by the reaction 
mixture. Hencx' as p is increased increases. 

Gelation envelope in the presence of difunctional monomer ; 

The broad prindplos discussed in' the previous sections 
will still hold in explaining the influence of the presence of 
B2 monesmor or when d 0 , Firstly it has already been noted 

that the vs, f' curve is not influenced by d and hence we 

need to discuss only vs. l . The relevant parameters 

therefore art' K and only. The limiting results for K and 
Kr— ^ O and CO can easily be derived and will not be presented. 



MT:r;j:,Vor^ tho oifocts of ;<. and K arc unalturod by the prtisonco 

R 

; 1' ^>2 ^i'ld hcnco nctod not bo discussed. Ww. shall focus our 
attention on the vffects of d. It is convenient to recall that 
' =. ( 1 -d) ( l-b) by ( 2B) , Firstly it should be noted that 
tiK.’ raaxii'num permisoiblo value for is ( 1 - d) when b = O. 
Further if d is increased while is kept constant, it implies 
a ds-crcasc in the concentration of the blowing agent in the 
reaction mixture. The results therefore can be understood in 
terras of the relative offoctivonass of B 2 and CD in creating 
connections between branch points. In the following discussicn 
it is assumed that j' is kept constant, 

( i) K = Oi Here CD is a nonreactive diluent. Thus if 

d'is increased, tlio fraction of B 2 monemor, which consumes A 

groups, is increased. We can therefore expect R . to increase, 

min 

This is shown in figure 7, 

Cii) K — > 00: Hero CD is an instantaneously reacting 
blowing agent and h^nco consumes more A groups than B 2 monomer. 
Hence* as d is increased, the slow reacting B 2 replaces the fast 
reacting CD resulting in a decreased This is shown in 

figure 8, 


( iii) = O; The CD monomer is a monofunctional blocking 

agent which consumes A groups wastefully. Thus as d is increased 
WG expect to decrease as B 2 does not have any blocking 

effect. This is shown in figure 9, 


( iv) Kj^-- -^oo: Here the CD 
in chain extension and providing 
the fact that the D group reacts 


mononujr is more effective 
linkages than monomer due to 

Z / 

at a faster rate than a B group. 
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FIGURE 7 : Effect of the addition of the difunctional 
rnonomer on Ks=0, Lines 1 and 2 

represent d « O and .5 respectively 
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PIGUllE 8 : Effect of addition of the difunctional 

monomer for K « oo . The curves 1-3 are 

for d = o, 0.2 and 0.5 respectively 






H.^nco v;o can expect to incr»jase as d is increased. This is 

shown in Fig, 10, 

For values of K and in between the above discussed 
linitS/ the behaviour is illustrated in figures 11 and 12 
and oan bo explained along lines similar to the ones used 
above and in the previous sections. 

Gel Point Calculation for Incomplete Conversions 

The curves fi.r incomplete conversions will lie inside 
the gelation envelope. If a value of C d 1) is fixed, all 

•k 'k ^ ^ 

thki other values o£ * 5 ^ calculatea 

from the kinetic parameters. For these probabilities, the 
values of R,; and d that will satisfy ( 35) and (36) can be 
dotcmiincd and plots of R versus r can be generated , It was 
found that the effect of rate constants on such plots was 
similar to that observed on gelation envoi opes, A typical 
plot is shown in figure 13, The curves start inside the lower 
half of the gelation envelope. To start with ^ 0*5 and it 
steadily, increases as R ie dtiangod. At scxne value of R, p attains 
a value of unity. This particular value of R will lie on 
vs,, ourve. ■ Such parametric curves give an indication of 
the extents of conversion attained when gelation occurs for any 
initial composition specified by R, and d lying inside the 
gelation envelope. 



FIGUi® 10 ; Effect of .addition of the difunctional 

monomer on for K.^ = oo * Tlie curves 

1-4 are for d *= O, 0,2, 0,5 and 0,8 
respectively. 




Effect of addition of the difunctional 
inonoitier on » 0*1 and = 1 

The curves 1-4 are tor d = O, 0,2, 0,5 
and 0,0 respectively. 
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FIGURE 13 » Gelation envelopes for inc3onplete conversion 
of all the monomers* The figure is for 
K w 0.1, K„ « 2, « 2,0, « 1.5 ana 

K« » 3.0, 

m 


CHAPTER 3 


In the previous chapter the reactions of formation 
of allophanate and biuret were assumed to be negligible* 

But at higher temperatures these reactions are quite 

2 €) 

competative ' , Even otherwise it is necessary to look at the 
effect of these reactions to determine whether their effect 
on gelation is really negligible or not. 


BASIC REACTIONS 


Isocyanates can react with the free hydrogen of urethan* 
and urea linkage to give allophanates and biurets respectively 
in the following x^ray: 


H 

! 

HO-R, 0=C~N 

1 ( ^ 1 . 
R^-C = O + R^-^=Cr=0 R - N - C 

a 



( 1 ) 


urethane isocyanate Allophanate 

H O H HO 1^0 H 

+ R2-W=CaO R-N-C-N-<^-lir-R2 C 2) 


urea isocyanate biuret 


It can be seen from the above equations that the allophanate 
reaction still leaves a free hydrogen for possible further 
reaction. But this kind of reactions do not take place in the 
industrial samples checked^, it is possible that the reaction 
stops at this stage due to steric hindrance. This also 
eliminates the possibility of the free hydrogen of biuret 
reacting further since these hydrogens are similar to the 
hydrogens of the allophanate group. 
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MODEL OF TOE SYSTEM 

A syfjtan consisting only of a dtisocyanate (A^A 2 )/ a 
polyol ( 3^) and water (CD) is oonsidored. The kinetics of 
each reaction is described below. 

Reaction of isocyanate with polyol and water has already 
been described in Chapter 2 and the same is adopted here. The 
syiT^bols for kinetic constants are also the same. 


Allophanate link is produced by the reaction of 
urethane with isocyanate. It is assumed that the rate of 
reaction differs from that with -OH by a proportionality 
constant A sample reaction step is as follows: 


A 4* B 


E ( urethane) 


E + A^ M ( allophanate) 


(4a0. 


E + ^ ^ allophanate) 


(4b) 


Let 0 be the fractional conversion to allophanate. It is 
defined to be the ratio of Allophanates formed to the total 
Allophanates that can be formed if all the B* s reacted 


(E)p - (E) 
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Blurot 

Tho biuret roaction rate constant is assumed to differ 
frora the urethane roaction rate constants by a proportionality 
constant of and the reaction scheme is as follows: 


+ D 

Kkj^ki 

F 

( urea) 

( 6) 

+ P 


G 

( biuret) 

(7a) 

+ P 

*f’^2 

G 

( biuret) 

(7b) 


Let be the fractional conversion tx.; Biuret, It is defined 
to be the ratio of Biurets formed to the total that can be 
formed if all the C's reacted. 

By writing the kinetic equations as before and solving 
thorn we can relate ^ and <0- toC<^ and ^ with 

reaction constants as parameters (See Appendix B) . 


staoichiomotry 

According to the above reaction schone: 

7v* s reacted = B* s reacted + C's reacted +D' s reacted +E's reacte 

+F' s react€= 




L.H.S. = A 


1 +oC + K*(l 


= B^3p + C^“*r+ C^?+ 


( 8 ) 

(9) 


Hence we obtain 
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R 


1 



+ ( i-^)(f+£+ <^) 

( 10 ) 


gelation criterion 


5 

Maoosko' s approach is used below to determine the 
average molecular weight of the polyurethane. This leads to the 
gelation criterion after finding the condition for. infinite 
molecular weight. 








( 11 ) 


As can bo seen above, B, E and E can lead to a branch, 

E(W = P( Abreacts) r P( B) react' ^ ^”b"' 

^^2 >- 

+ P(C) E(W^") + P(D) react'^^”^' 

+P(B)P(E^eaots’[®‘”“' + 

+P(D)P(P^,^„) \e(W^") t 

+P(E) E(wi") + P(F) E(W^’f^) 

+P(A2 not react) (O) 
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= i?( A , 

4'A( B) 

+?( E) 
Total B 


rcocts) 


r” 

'p(B)E(W^^) + P(G)E(«^^) + P(D)E(wff) 

3 C D 


T.PUt, 


4. + P(D) ?(F ^ 

reacts ' 2 ^ reacts F 


E(Wg^) + P(F) ECWp’^) I 

J 

tjquivalurit reacted = (f+i) + bJ£-+ C^0 


B 


E = “o'P 


B-(544) + C (f+ ^ + 5^) 


=(Bq + 2C, 




+«■) + 


Li 
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?(B) -ff^/ 

P(C) = LL-^ / Bj; etc. 
E(W«*^) =M„ + (f-1) E(WS^^) 


'B 


,in 


E(w 5“) = Me - + 


E 




D 


E(W^") - Me - M^^ + E 


’B 


( 12) 


( 13) 

(14) 

(15) 

(16) 

(17) 

(18) 


E(V^^^)= P(A,)E(W^) + P(A2)E(W^^) 

+?(A,)E(wJ;”) + P(A-)E(wJ”) 


A. 


2 J 


E(w5^^) = P (E 


reacts 


e(wS'^^) = P (F 


reacts 


-do- 


-do- 


1 


?( E 


) =e/p 


reacts 


P(F j. ) " 
'■reacts 


(19) j 

(20) I 

i 

f 

j 

(21) I 

i 

(22) I 

(23) ^ 



All 


- 1 * 

s roactod = reacted + A2 reacted + reacted 


+ A^ reacted 


A 


h ^ ( i^K, ) + i 

(i+K^oc,) ^ °(i+K,a:,) 


'A^l^ 

cK. 


• (“^2 + ■'a 


<^l) 


'■A''"! 


° (1 + 


1 4 0<2 + 1 


P( A ) tt „ ^ 

1 + CK. 2 t 1 t ^ 


( 24) 

( 25) 




P(,\*) « ^ ^ 


1 +o(^ + K^d + 0 <^) 


etc. 


( 26) 


E( 




" "■" "" " * ~~ — 

l+o <2 +^A^ ^ 


E(wj;’^) + Kp<J E(wj5? ) + K^ECW^^) 


'A* 


+C <2 E(W^" ) 



E(W^^'^) = 



X 


[ 

t 

€ 

P 


X 

X 


ECwJ’^)* E(W^) + :^(W^) 

tr i=> f > 

E(W^) « E(wJ^) + E(W^) 

ECW^’^) « PlA^) E(W^) + ^(A^) E(W^^”) + .,, 

1 2 


( 27) 

( 28) 
( 29) 

C 30) 

( 31) 

( 32) 

( 33) 
< 34) 


X 



Now using the ^ve aquations in 12, we obtain 


.^out 

/ % 




r* ) 





Mq+( f -1 ) 


f-fp 


U 


A 






(^■“c -«C02> + 


J, 


+ f6 X + Liril- J. X + (£-1) p X 




1(»C - +'(x 


+ X 


o<. 


A + B X 




( 35a) 


where 


A 




B 


+ 

C i-f) - ^ 
p {p +e) + -— (f+W i2f) 

(£-1)/^^ (jS+e + 2C^ + 2?+ g^^) + 2fe 

Li^ (f+i+ g^) 


j® (^+6-) T 2 

Substituting equation 35 in 27 we get 


( 35b) 


( 35g) 


P 


r 


f^K.W +-3^ ( A+BX)j + Q<^(m^ + A+BX') 


A<.‘ 



'k 

2A. + K o( ) 

M^+ 1 b-r 

^ I 4-0^^) 

2B (0<* + ) 

***** — --- — - ■ 

1 + K. (1 +o<.) 


( 36 ) 


Molecular woxcjht will be infinite if X is infinite hence 


gel criterion is 


k "k 

2B(o(o + K,X,) 

1 -fr '-'^2 4- 1 


>/ 1 


( 37) 


Of Chapter 


Without secondary reactions, £ = 0 = 


Hence, 


(f-D/’jg' 


l-f) 


(f+S) 


Hence (y'l criterion is; 


+ (!■/)& ic/^ + K^^) ^ ^ 

t ’'a' 1 +‘\'j 

2(f-l)M^ (Ol* + KfC^)^. (H.(<*)+ K^d +«) -2(1-/°) 

+(i=^)('( + ? 

4(f -l)/'^^ ( «■* + 

2[/A+ +‘^ + Vl+Vl- «1-/)S (<^ + K^'X^) 


Which is same as before hence the results are consistent. 









RESULTS AND DISCUSSION 


Using tho expression obtained for the gelation criterion 
the gelation limits were calculated and they are presented in 
Fig, 14 and 15, 

R < 1 

Effect of K^: As can be seen in Fig, 14 the gelation 

limit ^^Yitin as increased, Thxs happens 

b(2cause as Kg is increased more Allophanate links are formed 
and the number of crosslinks increase, causing the gelation to 
occur at a lower conversion. Hence the gel limit decreases both 
for a given R as /^^lin <3ecreaseS/ and for a givai as 
decreases. 

Effect of KgS, Increasing has a similar effect as 
increasing Kg at low values of , But at f = 1 since there are 
no urea linkages biurets are not formed and hence changing 
Kg has no effect. Also the effect of is very small for 
small K and values, 

R > 1 

As sQ&n in Fig,15 for nonzero Kg and values the 
limit for gelation increased for a given f . This happens 
because, the Allophanate and biuret linkages consune A groups 
and make gelation possible even at a higher A group concentration, 
in absence of such linkages such an excess of A (Above the 
lower line in Fig,15) can block the polyfunctional B groups and 
make gelation infeasible. 
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Thus it is clear that the presence of secondary 
rcacti;.'ns lead tt^ gelation of reactant mixtures which would 
not have been expected to gel. This is reasonable since 
these reactions ( at R < 1) provide more crosslinks for the 
same number of B groups while ( at R ^ 1) consuming more A' s 
for the same number of B's. It has been shown how even with 
very small reaction rates these reactions can affect the 
limits oonsiderably. Hence it is necessary that these 
reactions be taken into account while predicting the network 
properties. 



fi; 


It 1.8 known that ge? time measurement is a very good 
method for determining the consistency and uniform quality of 
the monomer batches used in industry. Consistency refers to 
the constancy of the ratios of monomers fixed for the particular 
system ^.;hilc uniform quality refers to the constancy of the 
purity of the chemical constituents of the system. Geltime method 
is applied to test these because geltime is sensitive to both 
reactant ratio as well as the reaction rates, 

17 

Typical observed gel time curves are shown in Fig. 16, 

The two curves v^ore obtained for two different systems 
as indicated in the figure. These figures show some interesting 
characteristics. Firstly there is a minimum geltime which may 
lie cither above or below a NCO/OH of 1. The second property 
of those curves is that they show a lower and an upper values 
of NCO/OH where gel time is infinite. These ratios will be 
referred to as lcw«r and upper gel limits. Beyond these limits 
gelling is not possible as gelation occurs when deficient 
reactant reacts completely. It is the purpose of this chapter to 
show that these two characteristic properties are explained by 
the nodels proposed in this work. The curves in Fig. 16 were 
obtained for solid elastomer systems where no foaming 
agent was used, 

MODEL FITTING 

First an effort is made to use an equal reactivity 
model for diisocyanates^ diol and triol to fit this data. The 
geltime curve for the equal reactivity case (with^ = 1# i.e. 



GSL T£MS ^ S<?c4.r!4s 



INDEX, NCO/OH 

' * * li it 

FIGURE |t6 I "Expeirlmental geLtJirie curves (Ref, 17), 

^ .‘iM Isocyanate used in FAP1@, The Polyols 
’ '“used are Ethyl trimetiiylol propane and an 
' ' ’ ' iaoilfOl,. 


Gel td»e 0 - 



NCO/bH 

FIGURE 17 1 . Calculated geltime curves. Is Equal 
reactivity case, 3s K,=si, K***!# 

K^j « 1.0, iCp « 0,1, 3; Kj^»l,K2*a2*5 

oY<^ina1& "no^" +0 5c«le . 


60 


no vjator) is shown in Pig, 17 (curve 1). The curve has the same 
overall s *rpo as experimentally obtained curves of Fig. 16, 

But a closer look reveals its limitations, While the 

position of the lower gel time is alright the position of the 
upper gel limit is at NCO/OH = 2 and the position of the minimum 
is at NCO/OH =1, Hence this model is inadequate for the 
genural system which may have a minimum at NCO/OH / 1 and an 
upper gel limit above 2, 

To ovojDcome the short comings of the above model, 
uncjgual reactivity of diisocyanate was introduced into 
the model. The results obtained are depicted in curve 2 o£ 
fig. 17, This shovJs a minimum gel time at NCO/OH y 1 and an 
upper gel limit at NCO/OH ) 2, This can be explained by the 
fact that infinite geltime indicates that at gel limits there is 
complete convorsion either of OH (R ')• 1, upper gel limit) or 
of NCO ( R ^ 1, lower gel limit) and as shown in Pig, 17, forc(^= 1 
and at « 1 ar upper gel limit of R ^ 2 is possible in 

unc 3 qual reactivity case, B’.’t *"his model failed to predict a 
minimum in gel time at NCO/OH ^ 1. inclusion of various 
possibilities (presence of moisture as an impurity or that 
1, functional group of Polyol / 3) did not alter this 
prc 3 diction. It was therefore found necessary to incorporate 
unaxjual reactivity of triols to explain this behaviour. Unequal 
reactivity of triols alone did not predict an upper gel point at 
NCO /oh ^ 2 as shown by curve 3, Fig, 17, Hence a model accounting 
for unequal reactivities of hydroxyl and isocyanate groups is 
required. The model is described below. 



DESCRIPTION OF MODEL 


The hydroxyl groups of diol are takai to be equally 
roactivu since tlic molecule is generally large and symmetric, 
the hydroxyl groups of the triol are assumed to have 
unequal reactivity of the structural asymmetry type. This 
triol can thus be represented as, 

% — ® ( 1 ) 

The reactivities o£ these and B 2 groups are assumed to diffe 
from the B groups (diol) by proportionality constants K and 
respectively* 

Thus# 


Aj + 

B 

> 

product 

( 2a) 

jAi *# Hh 

B, 

K 

product 

product 

( 2b) 

1 

1 

> 

( 2c) 

The method 

developed 

7 

by Miller and Macosko 

(also 


see Chapter 2) vjas used to derive the gel criterion. The 
following definitions were useds 


m 

0 

ll 

initial 

cone €21 tra toon of 

triol molecules 

( 3a) 

u 

0 

initial 

concentration of 

groups = 2B» 

{ 3b) 

&210 = 

initial 

concentration of 

B 2 groups = B^ 

(3c) 



initial csoncentration of diols 


B, 


( 3 d) 




2B2 = initial ooncoatration of B groups ( 3 e) 

o 


The concentration of A groups are represented as in Chapter 1 


R 




2A^/(2B2 + 3B3 ) 

o o 

3B3 /( 2B2 + 38 - ) 

o O ■^o 


( 3 f) 
( 3 g) 


The fractional conversion of B, Bj^ and 83 are defined as; 


( 4 b) 
( 4 c) 


5 = C 2B2 - [b] ) / 2B2 ( 4 a) 

/ o 

A =(283 

o o 

y^2 = ( 2B3 -[b^ ) / 2B3 

where [.®J# ooncentrations of B, B^^ and 83 

respectively in the solution. 

Using the same ^proach as that used in Chapters 1 and 
2, we obtain 


(1 = (1 -p 

U -f> 

and gel criterion is: 


K, 


( 5 a) 

( 5 b) 


1 ^ 2 y 


(gX* + 


( 1 +0(2 + 1 + 0 i ^) 


( 6 ) 


where 


Y = 


(1 -/)/> + *p2> 



ThG definitions of conversions of Isocyanate groups and their 
intorrelacions remain the same as in Chapter 2. The 
stt'ichium .itric relation thoi is given by 


fli 


1 +Q(* + K^(l -4^)^ 


2 ( 1 + 



(7) 


Thus the conversions at gel point and in particular can 
be computed using the above formulae. 


CALCULATION OF GELTIME 


Gel time at a givaiNCO/OH ratio was calculated by 
solving the first order differential equation 

- = k^ci + K^) (8a) 

which is the equation for consumption of diisocyanate monomer. 
This can be simplified in terms of the conversions as 


a|< = (i-f) + K{i -Pp 2^ + Kj^d -ji) 

! ( 8b) 

where, 

© = dimensionless time 

- -b£b( + 2B2 ) 


This equation with equation 5, 7 aid equations A13,14 and 15 from 
Appendix A was integrated by Runge-Kutta-Vemier method 
(IMSL subroutine) to a value of ^ 2 , satisfies the gel 

criterion 6, For a given sot of reaction rate constant 
ratios, geltime was determined at various values of initial 


) 



NCO/OH ratios and the position of minimum in geltime was 
graphically located* An example each where geltime minimiwi 
occurs at NOO/OH ^ 1 and ^ 1 respectively is given below: 


case 1 : Gel time bedow NGO/OH = 1 


Kj = 0,5, = 0*75, K2 = 3.0, K = 1,0, = 0,1, 

= 0,2 


NCO/OH 


0.5 

inf inite 

0.7 

7.16 

0.85 

4,74 

0.95 

4.56 

1.00 

4.59 

1.05 

4.70 

1.10 

4.86 

the minimum gel time lies below an NCO/OH 0 

Gel time above NCO/OH = 

1 

= ,5, = .5, = 

1.0, = 1 

NCO/OH Gel time { dimensionless) 

0.5 

inf inite 

0.6 

4,1 

0.8 

2.9 

1.0 

2.45 

1.5 

2.3 

1.9 

3.3 

2.0 

inf inite 
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To illustratG the position of the gel limits the 
following example will suffice: 

for = 4,0, K 2 = 4,0 and == 1,0, The position 
of upper gel limit is at 4,5 as seen frcim Fig, 6. 

It is concluded that inclusion of unequal reactivities 
of OH as well as NCO isnnecessary and sufficient to enablo the 
prediction of the two special characteristics of Fig. 16, experi- 
mentally known for polyurethane ejftstomers. 


/ 



CHAPTER 5 


In many polyurethane foam systems water is used as a 
blowing agent tB produce CO 2 on reaction with NCO, At low 
molecular weight of the mixture the viscosity of Polyurethane 
is low and CO 2 bubbles evolved will escape. At very high 
molecular weight the polyurethane is elastic and CO 2 evolved 
may elastically stretch the matrix whic>^ may contract after 
cooling, resulting in a change in dimension of the product. 

Hence there is a molecular weight range between which the CO 2 
evolved is effectively used. A discussion of this mechanism 
can be obtained elsewhere , 

In this chapter the intention is to calculate the CO 2 
evolved versus molecular weight. Such information is important 
because it points out how much CO 2 will be effectively trapped. 

calculation of CO 2 EVOLUTION AND AVERAGE MOLECULAR WEIGHT 

CO 2 is evolved according to eguation 1,2 (Chapter 1, 

equation 2). The cumulative amount of CO 2 evolved can be 

found frcxn the conversion of C and the rate from rate of 

consumption of The molecular weight can be calculated from 

the conversion of the functional groups using the approach of 
5 

Macosco et al„ as siiown below. For brevity the detailed 
steps are not repccrduced here. The formula is given below; 

= E(A) E(N^) + E(B^) E(Wg ) + £( 62 ) ECWg ) + E(C) E(W^) 

f 2 

( 1 ) 

( 1 + K + K OL) 

E(W,) = ^ ^ y 

^ a + ® 

E(Mg ) = f + Mg 

£ f 
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E(Wg ) =2 E(Wg^^) + M 




E( A) 
L 


ECB^) 

ECBg) 


E(C) 


ECWg ) 


B 


B. 


(l-T) +fj( - Meo^) + (!=► -) 


= R/2 L 


= R/2 + i/£ + d/( l«d) 

= ^fL 

ha 




2(1 - d) 


( 1 *" d) 2L 

i / -1 - d - X A * * 

= 2 ^ 1 - cl (^^2 + 

BM. +2AXr:: 


A il + + K^l) 

1 - 2^Y ^ ^ 


+ +CX 2 + ^ 


+ i i - =a - ril cf+i') 

(i-aXf-i) + dj + Li . ': - g - - . £ l(2^) 


X 


m 


- l - r < 3 -r, S . 


<A*^, Y , $ Gtc, are calculated by the equations derived in 
Appendix A, 

^ _ (1-C)^)/(1 + ^A*^l^ calculated from the following 
kinetic expression derived from the conservation equation of 
AA monomer (Equation A1 of Appendix A), 
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( 1 + K^) __ 

n '.'UY"^ 




r 


) + 


K(l-d - ) 


-f) 


^ (1-t) + 


( 2 ) 


where^ 

Q> ss t ( Bq + 2C^) = dimaisionless time 


is calculated from stoichiometric equation 1,36, 

Equation 2 was solved with Runge-^utta method to obtain 
conversions at different times and molecular weights as well 
a.a cx ) 2 evolved. 


discussion 

In Fig, 18 the CO 2 evolution rate is plotted against 
log (M^) for constant with R as parameter. From the figure 
it can be seen that as R increases the CO 2 evolved for a given 
increases. This can be easily ejqplained. A higher value of I 
implies a higher number of NCO groups per OH groups and hence to 
achieve a given i,e, same average number of branch points 
per molecular, a higher conversion of OH groups is required. 

This also implies that a higher number of water molecules have 
reacted leading to evolution of larger amount of CO 2 , 

In Fig, 19 R is constant while is the parameter for 
CO 2 evolution vs, log (^) curve. It is observed that as 
increases CO 2 evolved for a given decreases. This happens 
because bo achieve a given a system with higher requires 
a lower conversion. At the same time the amount of water 
available is low. Both these factors lead to a lower CO 2 
evolution for higher ^ value. 




HCRIWB 1)8 1 ®V61v6k 3 log a 

«• parameter. (/^« 0,5, K « 0,1). 

<i-/°; '^■/Z 



2 


3 


4 


( M^) 

FIC5URB 19 s CO 2 evolved versus log 

/* '"1 parameter. ( Rs=l,0,Ks5.0) . 
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USES 

Once sucli plots have been obtained they can be 
used to design ti*e system. For example in Fig*18 if = 1000 
is the most desiiable molecular weight at which CO 2 should 
evolve then R - 2,0 is the best value to use whereas R = 1.0 
or 100 is uns-itlsfactory™ “lb be able to use R = 1^0 we have 
to change ttie reactivity ratios in the appropriate direction 
using catalysts. 





coiclusiqns 

In Chapter 2 a model for describing the curing 
behaviour of a monomer with unequal reactivity effects has 
been developed and applied to systems consisting of a 
diisocyanate a polyol ( , a diol ( B 2 ) and foaming 

agent water (CD), The gelation condition was derived and it 
predicts { i) whether a given initial composition of the 
reaction mixture will lead to a gel or not and (ii) the 
fractional conversions of functional groups required before 
gelation can occur. 

Gel formation at complete conversion of the groups in 
stoichiometric deficiency defines the boundaries of gelation 
envelope and only those initial compositions which lie inside the 
envelope can lead to the formation of infinite networks. 

The boundaries of the gelation envelope have been presented as 
a plot of the ratio of the concentration of A groups and the 
concentrations of groups with which A can react, 

^in' the fraction of the polyfunctional groups f* • 

Water can influence gelaticai behaviour in three ways by 
acting like a non-reactive diluent, a monofunctional blocking 
ageit and an effective chain exteider linking branch points. 

As the water content in the reaction mixture is decreased, 
would increase if it were to operate like a nonreactive diluent 
and would decrease if the other two ways were operative. 

The total behaviour is a combination of all the three aspects 
and is determined by the relative speeds of reaction of water 
and hydroxyl groups. 



asymmetric characteristics of the monomea: 
can either promote or retard chain extension. Structural 
asymmetry retards chain extension. A first shell substitution 
effect that makes polymeric A ends to react faster than 
monomeric A ends promotes chain extension.. At a constant P ^ 
^ax found to increase if the overall effect of structural 
asymmetry and first shell substitution effect was to promote 
chain extension* 

Calculation of conversions at which infinite networks 
form for compositions which lie inside the gelation envelope 
has also been illusjxated. 

The calculations show (see Fig, 6) that an unequally 
reactive system may gel when an equal reactivity model 
predicts that gelation will not occxir (curve 4 for example) 
it may also not gel when equal reactivity model predicts that 
gelation is possible ( curve 2) , Hence it is important to 
consider unequal reactivity while making such predictions. 

The effect of secondary reactions has been shown to be 
not negligible (See Fig, 14 and 15) for R 1, presence of 
secondary reactions increases the limit fo^^ gelation. 

This means thati v;hile a model which does not consider 
secondary reaction may predict that system will not gel, in 
reality it may gel. For R ^ 1, secondary reactions lower the 
limiting means that the system which might be 

predicted to not gol may gel. Hence it is necessary to include 
secondary reactions in the model for correct prediction of 
gel limits. 
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I*- ’»jas shown that it was necessary to incorporate 
not only the unequal reactivity of NCO group but also that 
of triols to account for tl^ interesting observations made on 
gel time vs. N20/OH ratio, it was shown that asR increases CO 2 
evolved increases for a givai Also, as increases the 

00 2 evolved decreases. 



rtPI'ENDIX A 


If it is acsumed that reactivity of an end group is unaffe*- 
ctvti t:‘j dir.in length and that second order kinetics apply, 
tlifc following equations can be writt^. 


d ■A'J 

dF 

d aJ 
"'"c3t 


d A 


2 

HF 


{k^+k2) \A\ ; [b': + K [ c j + K f 


k, Yb + K’Ci+ K Kq i'd] ; 

^ I R - • - 


d!A 


i- - 


dA: 


dt 
d |>* J 


+ k^ rA/:'>’B]+ KiC')+ K K,, '"Dj > 


dt 

<3lCJ 

dt 


y * * ( . 

“4- k A • ’ I 
dt ^ 1 .1.' ,• 

dA. 


★ „ ':le 


K [c]U^1+5'4) + ^*2 ^ 


(Al) 

(A2) 

(A3) 

(A4) 

CAS) 

(A6) 




dl^Bj 
cb 


★ .. ^ 'k 


=1^ = 'B-;\(k,4.k2) tA]+ k^ ■.A3_> k^ fA^] V (*ve) 

r r . . , , , Ir 

since initially only nionaner is present, Aj|_ s= A^ = O initially 
Dividing (A3) by (A2) and integrating with the above initial 


condition 



K 


A 



(A9) 


where Combining this with (16) and (17) it is easy 

to show that 


~ ^2^ ^ 
o<.^C 1 

k k 

Initially 

f or ( A4 ) and ( AS ) ^5; 




= K, 


(AlO) 


-0^2^ 

= O and these are the initial conditions 
Suubstituting (A3) into (A4) and dividing 



the resulting equation by (Al), 


cl ' A, 


■ A- 


} 




lCj+k2 .A.5 


The initial ’'^‘^ticroor concentration is i Hence the above 

can be integrated^^witii the condition that fA^j = 0 when 

'A'-' ; 


A;^. The result is 


* 

[.Hi 


K, 


‘A] 


'O 


K., . -,-1 

A • ^ A 

K^-I-K^ ‘ ;A-t 
* 


- i ) A 

O V — O 


(All) 


where similarly substituting (A2') into (A5), dividing 

by CaI) and integrating 


* 

'.^2- 1 

.^ o ° ■.. l^-'o 


% 

lA / \ £4~1 

- C — ) 


■Ai 


( A12) 


Where K 2 = ^2^2* ’^^'^se results arc identical to those given 

in reference and hence aquations (20), (23) and (24) of 
that reference, ropredueso below, are also valid 


.A3 






A • 

1. -o 


1 ~'X, 




K^/l+1 


V 1 -X, 


2 *t H 


K. 




1 - Ai(K2-l-i^ / K^.Vl *■ 'V<1+1 


- (^^) ( 


1 -y2 = 




X 2( ^ l»K^) i ^a\+^ 


- ( 


1 -i-'> 


V2 




:!_-) K^t 1 


(A13) 


(A14 ) 


(A15) 


The above can be easily obtained by combining (16) to (20) 
with (A9) to (A12), For the other groups, the initial conditions 
are [Cj = O and = = 2 +rB 2 lQ. 

Dividing (A6) 1^ (A8) and integrating after using (21) and (22) 



t i-Kv 


(A15) 


Dividing (A7) by (A6) and integrating after using (22) and 
C 23) 

.. = •' . i ^ ' 1-(1-/) ^ (A17) 

1 L. ^ — 

Substituting these relations in ( 24) and equating it to ( 30) 
it is easy to show that stcichiometric relation reduces to 


O'.I 


iff if 

\ 1 + + K^( 1 4C-^) 


K^^^l+l ; ^ ^ "2 ‘'A' - -1^ '= R( 1 - d) i 

which is ( 36) • 


! (1 ^)C S + ) 


Similarly it is easy to show that 

if it 

v,> 

p. 


BB 


V.X 

= 4 5. 


2 


1+X+ 1 




2-,t:+(l -d-.«)(Y + S* ) 


CD 


= 4 


1 


2 1 
( 1 + 1 + 


( 1 -. d ~ r) ^ 


2j^|g+ ( 1-d-O i'i+i ) 


-' ( 1 - d) 


,/3 


2f'A4-(l - d -P) ( 'I +^) 


These relationships can be substituted into the gelation 
condition (34) tc obtain (35). Further if any one probability 
is known the rest can be calculated for a given initial composi- 
tion specified by R, pand d. For example if :V(i known ^ ck 

it 

and r /.^2 calculated by (A13)^ CA14) and (A15) respectively* 

/ c 

and /j are related to 0- and hence the stoichiometric relation— 
ship ( 36) can be used to obtain , The value of at any 

12 

time can be obtained by integration of (A2) as indicated earlier , 
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APPENDIX B 


The complete set of reactions can he represented by 
the following equationst 

- = ( V'2’rate+ kLc1+ k Kj,1pJ+ KgLE]+ Kpljlj 

(B.l) 

dfAfg 


at 

dCA^J 

dt 

dtA^I 

at 

dlA*l 

at 

d[cl 

dt 

dM 

dt 

dLBj 

dt 

d Ce3 

dt 


- LAj 


= k. CaI 




V 


dCAjl 


)7 


2 


+ k 


1 \ 


dt ■" “1 '- x-^ L 

dt 


V 


P . t; tA;i ^ ” 

kg[A23|K|c3 


} 


= S 


i 




;; 


j* K Kj^ D 

5 LBJ 


1 


J> - KgLEj) 


3^ = J)l{ -KKrKI + KptF] ) 


Using egns* Bl-BO we obtain 

,K 


(B.2) 

(B.3) 

(B.4) 

(B.5) 

( B, 6) 
(B.7) 

(B.8) 

( B.9) 

(B.IO) 


1-f = (1 -f3)' 

( = f- 

The equations relating and 0^2 iare sam€^ as A14 and A15 of 


(Bai) 

(B„12) 


Appendix A 



From equations B9 and B8 we obtain 


d 


1 


= 1 


K 


w 


ith boundary condition "[jeJ = ^ when ][b3 = 


we obtain 
o 


w/ 


1- LV 


iBl V 


1 K, 


Iji -p 


1-K, 


-1 


( B.13) 


Or, 


(3 -6- i„^Kg "P 


Or, 


= r ( 



1 


From equations BIO, B6 and B7 we obtain? 


a W 

d (C)‘ 

d[D] 


d IdI Cli /V 

Hh T-/*y •n ^ ry ' 


<a icj XC2 


■K 


K 


R 


Kr -1 


K. 


1 . (isL, 

_ L'^ol 


_fs_ fi _ 

-1 I 


K 


(B.14) 


d (C) 

for consumption of D only. 

Or, 

M _ IfI , ‘S' , 

-TtCJ - (CJ K 
solving with boundary conditions [f] = O at |cj = we 
obtain 


R*^l 


c^y Co 


K 


R 


‘"o] 


in'fr r-m i-,'Twmt-i I i ittm 

■% “ ‘'r "1 


1 1, 


k-r-Kj 1 
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- Kj, 




Or, 


9 ^ 


C . 

'(1 4 

® -1 

R 



(1 »f) (1 - ( l4) 


R 




K 


R 


(B,15) 
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APPENDIX C 


Macosko’s method was used to derive the equation 
for sol fraction and network density. The formulae derived 
are reproduced below. For brevity their derivation is not 
presented here 


Sol-fraction: 


For Model 1 ( as in Chapter 2, 1) 


= sol fraction = Dl -i- d2 + d 3 PAl PA2L 

o 


+ D4 PA2 PAIL + D5 


(Cl) 


ok = (1 + 


( irjdo 1 
^ d “ d 


1 4 2 ^ ( 1 ^ ) 

■ 0(*2 ^ + 1 + I <;^ 


B 


1 + a<2 + ( 1 

Ll + +0^2 + 


/d 


c = ( -A 4 - N - 4AB) /2 


dC^ + (1 -^)C 


D = 

PAl = D 
PA2 = D ' 

PA2L = ( 1 +««^D = pail 

DO 

Dl 
D2 


= d /3 + ( 1-1) M /2 4- RM, /2 
3 2 2 

= a Mg /3D0 


= (i-f) 


Mg / 2D0 



D3 

D4 

D5 


^ ^ y.) /do 


1 .;. 


K H, /DO 


(•— ■)(<X) H /DO 
2 


Avarago i-iol.icular woi./at between branch points (M^) 


For Model 1 (Chapter 2, No water) 

P, 


f. 


B 


(1-Pnn^ ^ ^ 


P P 


BB B£ 


2 ( 1 “ 


] 


B, 


BB 


2 ( 1-^) 


o<; + 


1 

L 1 +<?<^ 


2 + ( 1 


r c^2 + 


2dp I -- 


i ^1 i-O^* + K , (1 +<X*) 


(C2) 


Network density was derived using Macosko’s approach for 
Model 2 (Chapter 3) xnr the case of no water present. 


Network concentration 


+ 3c-(i-ip (1 - 


(3Fb-i+^) + ^1^(1 - Fb>(1 - *'a> ^b‘*'b - 1 +f> 
+3C1 - (Fg - lY)^eVfi3 ^ 


Fq = ( 1 -'■^1 - 4 MN) /2I'1 


= B + (p2^ + j5C-) C/(|i+6(l - O) 


p + pc-) c/(p^ 


A 

M = pc ( p+e) / ( p+eci “ c)) 

N = (1 -p) + b (g+<r(l - C)) 

(1 +1^) - (c4 + 


B 


1 +0<x + K, (1 +0<, ) 


C = 2 


2 

cx; t 


1 +C<2 + ( 1 + 


(C3) 
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